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Abstract 

The production of photons accompanied by jets in large-Q 2 deep inelastic ep scatter- 
ing is calculated in next-to-leading order. We describe how to consistently include 
contributions from quark-to-photon fragmentation and give numerical results rele- 
vant for HERA experiments. 



1 Introduction 

The production of isolated photons in high-energy hadronic collisions is an important 
testing ground for QCD. Since the photon does not take part in the strong interaction, it is 
a "direct" probe of the hard scattering process and provides a means to measure the strong 
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05 7 HH 92P (0), and by EU Fourth Framework Program Training and Mobility of Researchers through 
Network Quantum Chromodynamics and Deep Structure of Elementary Particles under Contract FMRX- 
CT98-0194 (DG12 MIHT). 
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coupling constant a s or to extract information on the parton distributions, in particular 
the gluon density in the proton [ij. Moreover, good knowledge of the standard model 
predictions for direct photon production is required since it is an important background 
for many searches of new physics. 

At HERA, with increasing luminosity, the measurement of isolated photon production 
will give information on the parton content of the proton and at Q 2 = 0, i.e. for pho- 
toproduction, also on the parton distributions in the photon. First experimental results 
from the ZEUS collaboration 0, |3j (see also H) have been reported and found in rea- 
sonable agreement with next-to- leading order (NLO) predictions |J. Cross sections for 
the production of hard photons in deep inelastic scattering are much smaller as compared 
to photoproduction and, therefore, more difficult to measure. Typical cross sections, for 
example with Q 2 > 10 GeV 2 , are of the order of 10 pb. With a luminosity of 50 pb -1 one 
thus expects measurements of differential cross sections to become feasible. 

A NLO calculation for direct photon production in deep inelastic ep scattering, ep — > ejX, 
at large Q 2 has been reported recently by two of us and D. Michelsen ||. Since hard 
photon production occurs, compared to inclusive deep inelastic scattering, at a relative 
order 0(a) one expects a sizable cross section only at moderately large Q 2 . Therefore, one 
can restrict the calculation to pure virtual photon exchange with Z exchange neglected. 
In H the hadronic final state is separated into 7 + (1 + 1)- and 7 + (2 + l)-jet topologies 
(the remnant jet being counted as "+1" jet us usual). The approach is thus analogous 
to the calculation of (2 + 1)- and (3 + l)-jet cross sections, where one of the final state 
gluons is replaced by a photon 0. In addition to the direct production, photons can also 
be produced through the fragmentation of a hadronic jet into a single photon carrying 
a large fraction of the jet energy ||. This long-distance process is described in terms of 
the quark-to-photon and gluon-to-photon fragmentation functions which absorb collinear 
singularities present in the perturbative calculation. First measurements of the q — > 7 
fragmentation function in e + e~ — > 7 + 1-jet are presented in || (see also []TD| for the 



discussion of an inclusive measurement). The NLO theory for this process has been 



worked out in JTTJ . In |J the fragmentation contributions were discarded and the photon- 
quark collinear singularities had been removed by explicit parton-level cutoffs. The results 
depended strongly on these photon-parton cutoffs, in particular for the incoming gluon 
contributions [[J. However, these cutoffs are difficult to control experimentally. 

In this paper we report results in which the fragmentation contributions are included 
together with isolation criteria that limit the hadronic energy in the jet containing the 
photon. Whereas in || various photon + jet cross sections with invariant mass jet reso- 
lution criteria were calculated, in this work we concentrate on the calculation of various 
differential cross sections either exclusively or inclusively which depend on the transverse 
momenta and rapidities of the photon or the accompanying jet. We use the 7*p center- 
of-mass system to define the kinematic variables. The cone method is applied to define 
the parton jets and to isolate the photon signal. 

In section 2 a brief outline of the theoretical background for calculating the cross section 
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is given. In section 3 numerical results are presented. Section 4 contains a short summary 
and the conclusions. 



2 Subprocesses Through Next-to-Leading Order 
2.1 Leading-Order Contributions 

In leading order, the production of photons in deep inelastic electron (positron) scattering 
is described by the quark (antiquark) subprocess 



where the particle momenta are given in parentheses. The momentum of the incoming 
quark is a fraction £ of the proton momentum pp, p% = ^pp. The proton remnant r has 
the momentum p r = (1 — £,)pp. It hadronizes into the remnant jet so that the process (1) 
gives rise to 7 + (1 + l)-jet final states. In the virtual photon 7*-proton center-of-mass 
system the hard photon recoils against the hard jet back-to-back. To remove photon 
production by incoming photons 7* with small virtuality (photoproduction channel) and 
to restrict to the case where the scattered electron e(p2) is observed, one applies cuts on 
the usual deep inelastic scattering variables x, y and Q 2 . In addition, to have photons 
"f(ps) of sufficient energy we require an explicit cut on the invariant mass W of the final 
state, W 2 = (q + pp) 2 , where q is the electron momentum transfer, q = pi — p 2 and 
Q 2 = —q 2 as usual. Both leptons and quarks emit photons. The subset of diagrams 
where the photon is emitted from the initial or final state lepton (leptonic radiation) is 
explicitly gauge invariant and can be considered separately. Similarly, the diagrams with 
a photon emitted from quark lines is called quarkonic radiation. In addition, there are also 
contributions from the interference of these two. Since we are interested in testing QCD 
under the circumstances that the photon is emitted from quarks the contributions from 
leptonic radiation are viewed as a background and must be suppressed. This can easily 
be done by a cut on the photon emission angle with respect to the incoming electron ||. 
In our numerical evaluation we include this background source as well as the interference 
contribution. 

At lowest order, each parton is identified with a jet and the photon is automatically iso- 
lated from the quark jet by requiring a non-zero transverse momentum of the photon or jet 
in the 7*-proton center-of-mass frame. Therefore the photon fragmentation contribution 
is absent at this order. 

2.2 Next- To-Leading Order Corrections 

At NLO, processes with an additional gluon, either in the final state or in the initial state, 
must be taken into account, i.e. 




(1) 
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e(pi) + q(p 3 ) -> e(p 2 ) + g(p4) + t(Pb) + 0(Pe), (2) 
e(pi) + 0(pa) -> e(p2) + g(p 4 ) + 7(Ps) + g(pe), (3) 

where the momenta of the particles are again given in parentheses. In addition, virtual 
corrections (one-loop diagrams at 0(a s )) to the LO process (1) have to be included. 



The complete matrix elements for (2) and (3) are given in [12]. The processes (2) and (3) 
contribute both to the 7+(l+l)-jets cross section, as well as to the cross section for 7+(2+ 
l)-jets. In the latter case each parton in the final state of (2) and (3) builds a jet on its own, 
whereas for 7 + (1 + l)-jets a pair of final state partons is experimentally unresolved. The 
exact criteria for combining two partons into one jet will be introduced when we present 
our results. Following the customary experimental procedure the resolution constraints 
will be based on the cone algorithm. 

In the calculation of the cross section for 7 + (1 + l)-jets we encounter the well-known 
infrared and collinear singularities. They appear for the processes (2) and (3) in those 
phase space regions where two partons are degenerate to one parton, i.e. when one of the 
partons becomes soft or two partons become collinear to each other. The singularities 
are assigned either to the initial state (ISR) or to the final state (FSR). Contributions 
involving the product of an ISR and a FSR factor are separated by partial fractioning. 
The FSR singularities cancel against singularities from the virtual corrections to the LO 
process (1). For the ISR singularities, this cancellation is incomplete and the remaining 
singular contributions have to be factorized and absorbed into the renormalized parton 
distribution functions (PDF's) of the proton. 

To accomplish this procedure, the singularities are isolated in an analytic calculation with 
the help of dimensional regularization. This is difficult for the complete cross sections of 



the processes (2) and (3). After partial fractioning, the phase-space slicing method |13 
is used to separate the singular regions in the 4-particle phase space and to determine in 
these regions the approximated matrix elements and phase space factors. In those regions 
only, the calculation is performed analytically. For this purpose a slicing cut y$ is applied 
to the scaled invariant masses yij, where y^ = (pi+Pj) 2 /W 2 with W 2 = {pp+pi—p2—ps) 2 - 
j/o must be chosen small enough, so that terms of the order 0(jJq), which are discarded 
due to the singular approximation, are so small that an accuracy of a few percent can be 
achieved for the final result. 

To be somewhat more explicit, let us assume that by partial fractioning the contribution 
proportional to the pole term I/2/46 has been isolated in the matrix element \M\ 2 for the 
process (2). In the infrared region p§ — > the two invariants y^ and y$$ vanish. Then the 
integration over these two variables is performed: (z) over the singular region (S), y 46 < y$ , 
?/36 > 0, analytically with 4 — 2e dimensions and using the singular approximation; (ii) 
over the finite region (F), y m > y 36 < y£, numerically without any approximation 
to \M\ 2 and in 4 dimensions; and (Hi) over the explicit two-parton region (R), y iQ > y$, 
1/36 > Voi s ^ so numerically. This separation yields two contributions, the parton-level 
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7 + (1 + l)-jets, which come from the integration over the regions (S) and (F) and the 
parton-level 7 + (2 + l)-jets contribution, which corresponds to the integration over the 
region (R). All the remaining phase space integrations are performed numerically with 
usual Monte Carlo routines. For the physical cross sections defined in the next section, 
which are obtained by adding the contributions from the regions S, F and R, of course, and 
after adding the virtual contributions and performing the subtraction of the remaining 
ISR collinear singularities, the dependence on the slicing parameter y$ cancels. This has 
been checked in 0. This means that the cut-off is purely technical. 

In addition, the matrix elements |M| 2 for the processes (2) and (3) have photonic infrared 
and collinear singularities, i.e. due to soft or collinear photons. The infrared singularity 
can not occur since we require a sufficiently large photon energy E 1 = \p<j\. But collinear 
singularities are present in general. In the earlier work || these g-7 collinear contributions 
were eliminated by an isolation cut on the photon of the form y 5i > yl with a sufficiently 
large isolation parameter yl which was considered as a physical cut. In this approach, 
the photon was considered as a special parton, which was always isolated from all other 
partons in the initial and final state. Such a photon isolation is very difficult to impose 
experimentally, since it refers to a separation of the photon from partons, whereas in the 
experiment only hadrons are measured directly which are recombined to jets. In addition, 
it was found in M that the results, in particular for the gluon-initiated process (3), depend 
strongly on the isolation cut y^. In (3) the photon can become collinear to two final state 
partons, namely q and q, which explains the stronger dependence compared to the process 
(2), where there is only one quark (or antiquark) in the final state. Although under the 
kinematical conditions assumed in || the contribution of the process (3) was only some 
fraction of the total cross section for 7 + (1 + l)-jets, the y^ dependence of the final result 
was undesirable. In a more systematic treatment the y$ dependence can be avoided by 
adding contributions from the quark-to-photon fragmentation function (FF). In order to 
achieve this we include the contributions to \M\ 2 from (2) and (3), where y^i < yl with 
i = 4 in (2) and i = 4 and 6 in (3). This leads to collinear divergent contributions which 
are regulated by dimensional regularization. The divergent part is absorbed into the bare 
photon FF to yield the renormalized FF denoted by D q ^. The additional fragmentation 
contribution, which includes the contributions from the region y 5i < y$ (i = 3,4), has, for 
example, for the process (2) the following form 



M 
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M 



y*q->qg 



(4) 



where D q ^(z) in (4) is given by [[14 



D q ^{z) = D q ^{z^ 2 F ) + 




( 



P 97 0)ln 



z(\-z)yZW 2 

/4 



+ z 



) 



(5) 



P qi {z) is the LO quark-to-photon splitting function 




l + (l-z) 2 



(6) 



z 
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and e q is the electric quark charge. D q ^ 1 {z^ 2 F ) stands for the non-perturbative FF of 
the transition q — > 7 at the factorization scale fip. This function will be specified in 
the next section when we present our results. The second term in (5), if substituted in 
(4), is the finite part of the result of the integration over the collinear region y^ < y^. 
As will be explicitly shown in fl5fl , the dependence in (5) cancels the dependence of 
the numerically evaluated 7 + (1 + l)-jet cross section restricted to the region y 5i > yl, 
investigated in || . The variable z stands for the fraction of the photon energy in terms of 
the energy of the quark that emits the photon. Suppose the photon is emitted from the 
final state quark with 4-momentum p' A = P4+ p§, then z can be related to the invariants 
2/35 and y 34 , 

2/35 1/35 /-X 

z = = (7) 

1/34' 1/34 + 2/35 

The fragmentation contribution is proportional to the cross section for 7*g — > qg, which 
is of 0(a s ) and well known. It must be convoluted with the function in (5) as indicated in 
(4) to obtain the contribution to the cross section for 7*g — > qyy at 0(aa s ). Equivalent 
formulas are used to calculate the fragmentation contributions to the channel (3) and in 
the case where the quark in the initial and final state is replaced by an antiquark in (2). 



3 Results 

The results presented in this section are obtained for energies and kinematical cuts ap- 
propriate for the HERA experiments. The energies of the incoming electron (positron) 
and proton are E e = 27.5 GeV and E P = 820 GeV, respectively. The cuts on the usual 
DIS variables are 

Q 2 > 10 GeV 2 , W > 10 GeV, (8) 
10~ 4 < x < 0.5, 0.05 < y < 0.95. 

To eliminate the background from lepton radiation || we require 

90° < fl 7 < 173°, 6 je > 10° (9) 

where 9 1 is the emission angle of the photon measured with respect to the momentum of 
the incoming electron in the HERA laboratory frame. The cut on 9 ie = (e(p2), 7(^5)) 
suppresses leptonic radiation from the final-state electron. The PDF's of the proton are 
taken from |1| (MRST) and a s is calculated from the two-loop formula with the same 
A value (Ajjg(nf = 4) = 300 MeV) as used in the MRST parametrization of the proton 
PDF. The scale in a s and the fractorization scale are equal and fixed to Q 2 . 

We are interested in the differential two-particle inclusive cross section E^Ejda/d 3 p-yd 3 pj 
at NLO (up to 0(aa s )), where (E,p) represents the four-vector momentum of the 7 or jet. 
The 7 + (1 + 1) -jet cross section receives contributions from leading and next-to- leading 
order and the 7 + (2 + l)-jet cross section from leading order only. In the latter case, only 



6 



7 + 2-parton-level jets contribute, while each parton including the photon build a jet on 
their own. The evaluation of the 7 + (1 + l)-jet cross section is based on two separate 
contributions, a set of two-body contributions, i.e. 7 + 1 parton-level jet, and a set of 
three-body contributions, i.e. 7 + 2 parton-level jets. In this definition of parton-level 
jets the remnant jet " + 1" is not counted whereas the photon is considered also as a 
parton, like q, q and g. Each set of contributions is completely finite as all infrared and 
collinear singularities have been canceled or absorbed into the proton PDF or the quark- 
to-photon FF. Each contribution to the 7 + (1 + l)-jet cross section depends separately 
on the slicing parameter yl . The analytic contributions are valid only for very small yl . 
Separately, the two contributions have no physical meaning. For the contributions to the 
7 + (1 + l)-jet cross section coming from 7 + 2 parton-level jets, a slicing cut y$ for the 
photon is introduced. After adding the photon fragmentation contribution this sample 
becomes independent of y^, i.e. also y^ has the status of a technical cut like yl . In the 
7 + 1 parton-level jet event sample the photon is isolated from q and q by requiring a 
non-zero transverse momentum of the photon in the 7*p center-of-mass frame. The two 
technical cuts yl and yl serve only to distinguish the phase space regions, where the 
integrations have been done analytically with arbitrary dimensions from those where they 
have been done numerically in four dimensions. These two parameters must be chosen 
sufficiently small to justify the neglect of terms proportional to yl and yl, respectively 
in the analytical part of the calculations. We found that yl = 10~ 4 and yl = 1CT 5 is 
sufficient to fulfill this requirement. 

In order to comply with the jet definitions in future analyses of experimental measure- 
ments the partons and the photon in the 7 + 2 parton-level jet sample are recombined 
to 7 + (1 + l)-jets using the cone algorithm of the Snowmass convention |T7|]. In this 



recombination scheme the photon is treated like any other parton (so-called democratic 
algorithm). In the j*p center-of-mass frame, two partons i and j are combined into a jet 
J if they obey the cone constraint < R, where 

Ri,j = ^(rii-Vj) 2 + (<Pi-<Pj) 2 - (10) 

77 j and <pj are the rapidity and azimuthal angle of the recombined jet. These variables are 
obtained by taking the averages of the corresponding variables of the recombined partons 
% and j multiplied with their p T values. The px of jet J is the sum of Pt,% and ptj- We 
choose R — 1. In some cases, an ambiguity may occur, when two partons i and j qualify 
both as two individual jets i and j and as a combined jet J . In this case we count only the 
combined jet J to avoid double counting. The rapidity is always defined positive in the 
direction of the proton remnant momentum. The azimuthal angle is defined with respect 
to the scattering plane given by the momentum of the beam and the scattered electron. 
One of the recombined jets may be the photon jet. To qualify a jet as a photon jet we 
restrict the hadronic energy in this jet by requiring 

^7 = E^l = 1 — e had > z cut . (11) 

PT, 7 + PT,had 



7 



Pr,7 and Pr,had denote the transverse momenta of the photon and the parton producing 
hadrons in this jet, respectively. For our predictions we choose (-had — ^had = 

0.1 H- 

Eventually the parameter e° arf (or z cu t) must be chosen in accordance with the experimen- 
tal analysis. Our results depend also on the choice of the quark-to-photon fragmentation 
function. The factorization scale dependent quark-to-photon FF of 0(a) is taken from 



14}| . It is the sum of two contributions, the solution of the evolution equation at this order 
and an initial FF at some initial scale /xq- The initial FF function and initial scale have 
been fitted to the ALEPH 7 + 1-jet data [[|. The factorization scale dependent quark- 
to-photon fragmentation function also gives a good description of the inclusive photon 



distribution as measured by OPAL [10 



With these definitions it is clear that in NLO the final state may consist of two or three 
jets, where one jet is always a photon jet. The three-jet sample, equivalent to 7 + (2 + 1)- 
jets in the notation of the previous sections, consists of all 7 + (2 + 1) parton level jets, 
which do not fulfill the cone constraint (10). 




5 10 15 20 25 30 35 40 45 50 5 10 15 20 25 30 35 40 45 50 



(a) MGeV] (b) prlGeV] 

Figure 1: pt distributions of photon (a) and jet with largest px (b) for LO (dotted line), 
NLO 7 + (1 + l)-jets (dashed line) and NLO 7 + (1 + 1)- and 7 + (2 + l)-jets (full line). 



In Fig. 1 and 2 we show our results for the pt and 7] dependence of the cross sections, 
da/dpT and da/drj, concerning the photon and the jet with the largest pt- In each figure 
we have plotted three curves, (i) the cross section in LO (dotted curve), which has only 
7 + (1 + l)-jets and is independent of the jet defining parameters R and e° had1 (ii) the 
NLO cross section for 7 + (1 + l)-jets (dashed curve) and (iii) the sum of the NLO cross 
sections for 7 + (1 + l)-jets and 7 + (2 + l)-jets (full curve). Specifically, in Fig. la we 
present da/dpx^, the transverse momentum dependence of the three cross sections (i), 
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Figure 2: r\ distributions of photon (a) and jet with largest pt (b) with labeling of curves 
as in Fig. 1. 



(ii) and (iii) for p Ta > 5 GeV. All other variables, in particular r/j, and Pt,j, are 
integrated over the kinematically allowed ranges. We see that all three cross section have 
a similar shape. The sum of the 7 + (1 + 1)- and 7 + (2 + l)-jets cross section is only 
slightly larger than the 7 + (1 + l)-jets cross section. Both cross sections do not differ 
very much from the LO cross section indicating that the NLO corrections are not very 
large. Of course, this is a consequence of our choice for the cone radius R. In Fig. lb the 
plot of da/dp T j for the jet with the largest pt is shown. The qualitative behaviour of the 
three cross sections (i) - (iii) is similar as in Fig. la. For the rj distributions we integrate 
over pT n > 5 GeV and pt,j > 6 GeV. The choice of two different values of minimal pt's 
for the photon and the jet is necessary to avoid the otherwise present infrared sensitivity 
of the NLO predictions. This sensitivity is known from similar calculations of dijet cross 
sections in ep collisions [|18[] and must be avoided. The cross section da/drj^ is plotted in 
Fig. 2a, again for the three cases (i), (ii) and (iii). The shapes of the three curves are 
similar. Here we have integrated over the full kinematic range of the variable r]j. Figure 
2b contains the predictions for da/drjj, where rjj is the rapidity of the jet with the largest 
Pt- In Fig. 2b we observe that the full curve, which represents da/drjj for the sum of 
the two jet cross sections, is shifted somewhat more to rjj > compared to the NLO 
7 + (1 + l)-jets cross section. The LO cross section (dotted curve) peaks more in the 
backward direction than the other two. Compared to dcr/drj^, shown in Fig. 2a the rjj 
distribution for the jet peaks at somewhat smaller rjj. 

For completeness we also give the results for the various cross sections in Table 1, separated 
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Table 1: Contributions to the cross sections (in pb) for 7 + (1 + 1)- and 7 + (2 + l)-jet 
production (see text). 



KjOlllTl Dution 


7 + (,1 + ij-jets 


7 + \* I J- J _ J e ^ s 


LO 


8.561 ±0.022 




s 


-32.164 ±0.026 






24.843 ± 0.034 


0.6349 ±0.0008 


R(<?) 


3.717 ±0.005 


0.3239 ±0.0007 


F 


3.429 ± 0.022 




D 


-0.599 ±0.001 




sum 


7.787 ±0.053 


0.9588 ±0.0010 



into the contributions S, R (for incoming quarks or antiquarks and incoming gluons), F 
as described in the previous section and the fragmentation contribution denoted by D as 
defined in Eq. (5). About 12 % of the total NLO cross section is due to 7+ (2 + l)-jet final 
states. The 7 + (1 + l)-jets cross section is reduced by about 9 % by NLO corrections. 



4 Concluding Remarks 

We have presented a NLO calculation for the production of photons accompanied by jets 
in deep inelastic electron proton scattering, taking into account the contribution from 
quark-to-photon fragmentation. This improves a previous work which suffered from the 
presence of parton-level cutoff parameters. The present consistent treatment allows for a 
direct comparison of our theoretical predictions with experimental measurements without 
being sensitive to uncertainties from unphysical cutoff parameters. 

We expect that the measurement of photon plus jet production at HERA will contribute 
to testing perturbative QCD. Moreover, our results add another piece to the set of NLO 
predictions of the standard model needed in searches for new physics. The calculation 
covers the range of large Q 2 up to several 10 3 GeV 2 . At even larger momentum transfers 
additional contributions from Z exchange become as important as pure 7 exchange to 
which the present work was restricted. 
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